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NOTATION 

The  symbols  used  in  this  study  are  defined  where  they  first  appear.  For 
convenience,  frequently  used  symbols  are  summarized  below. 

The  bar  and  tilde  generally  denote  a  vector  and  a  matrix,  respectively,  whether 
as  underscores  or  overscores.  Braces  and  brackets  are  used,  respectively,  to  denote 
a  vector  and  a  matrix.  For  example,  a  column  vector  U  can  be  written 


—2 


ua  =  {ual 


with  its  subvector 


[B]  Matrix  relating  strain  vector  to  nodal  displacements 

[D]  Elasticity  matrix 


Young's  modulus  of  elasticity 


{F} 


A  load  vector 


F  ,F  ,F 
x  y  z 


Element  nodal  force  in  the  x,  y,  or  z  direction,  respectively; 
forces  are  positive  in  the  positive  direction  of  x,  y,  or  z  axis 


G 


E 

2(l+v) 


shear  modulus  of  elasticity 


H,  , H .  ,H 
ix  iy  iz 


Weighting  coefficient  corresponding  to  position  along  Gaussian 

quadrature  points  5.  ,  n.  ,  or  £ . 

ix  ly  iz 


v 


«*  i  -i  -i  it  ii  A-  Mfi‘r 


I 


Moment  of  inertia  of  a  transverse  section  of  a  beam 


i 

AAA 

[J] 

|J| 

[Ke]  or  [K]e 

NNPE 

P 

[SK] 

TEMP (I) ,TE(I) 

TEMPF(J) ,TF(J) 

TLOAD 

TREF 

t 

U 

u,  v 


V,  Vol 

X,Y,Z 

a 


Subscript  indicating  nodal  number,  or  active  index 

Vector  having  unit  value  in  direction  of  x,  y,  or  z  axis, 
respectively 

Jacobian  matrix  of  coordinate  transformation 
Jacobian  determinant 

Stiffness  matrix  of  an  element  e 

Function  of  curvilinear  coordinates,  taking  a  value  of  unity 
at  node  i  and  zero  at  all  other  nodes 

Number  of  nodes  per  element 

Applied  pressure  on  an  element  boundary 

Stiffness  matrix  of  entire  structure 

Temperature  at  node,  I 

Thermal  expansion  force  at  node  I 

Global  thermal  load  vector  (entire  structure) 

Ambient  temperature,  or  arbitrary  referenced  temperature 
Plate  thickness 

Vector  of  nodal  displacement  for  entire  structure 

Component  of  displacement  in  the  direction  of  the  x  or  y  axis, 
respectively,  displacements  are  positive  in  the  positive 
direction  of  coordinate  axes 

Component  of  displacement  at  node  i 

Volume  of  a  given  solid 
Global  system  of  rectangular  coordinates 
Coefficient  of  thermal  expansion 
Shearing  strain  components 
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ABSTRACT 


A  finite  element  code  to  account  for  thermal  expansion 
in  a  solid  was  developed  for  the  Independent  Research  and 
Independent  Exploratory  Development  project  "Tribology  of 
Sliding  Surface  Bearings."  The  program  is  based  on  a  two- 
dimensional  model  using  a  second,  or  higher  order  inter¬ 
polation  function  in  the  element  space  that  will  allow  a 
diverse  temperature  field,  such  as  a  steep,  nonlinear 
temperature  gradient,  to  be  prescribed  in  a  solid  body. 

The  computer  code  has  a  definite  advantage  over 
certain  finite  element  systems,  commercially  available. 

Many  accept  only  a  constant,  or  averaged  temperature  input 
into  their  element  space.  With  the  new  capabilities, 
complex  thermal  mechanical  responses  under  severe 
temperature  gradients  can  be  readily  analyzed.  For  instance, 
the  hot  spot  in  a  ship’s  landing  deck  due  to  the  concentrated 
heat  load,  such  as  those  generated  by  high  temperature  jet 
exhaust,  can  be  more  realistically  represented  by  the 
elements  of  current  development.  The  element  mesh  size  and 
the  input  data  set  are  more  manageable. 


ADMINISTRATIVE  INFORMATION 

This  report  covers  work  conducted  for  and  funded  by  the  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center’s  Independent  Research  and  Independent 
Exploratory  Development  project,  Tribology  of  Sliding  Surface  Bearings.  The  work 
was  accomplished  under  Program  Element  61152N,  Task  Area  ZR  000  01  01,  and  Work  Unit 
2832-121. 


BACKGROUND 

PRELIMINARIES 

The  description  of  the  behavior  of  solids  under  the  action  of  heat  and  mechan- 

1* 

ical  loads  has  been  established  in  the  field  of  thermoelasticity,  however,  the 
determination  of  the  elastic  or  inelastic  deformations  in  solid  bodies  under  pre¬ 
scribed  temperature  distributions  is  only  feasible  for  a  number  of  elementary  cases. 
It  is  difficult  to  generalize  to  more  complex  configurations  commonly  existing  in 
real  structures.  Practical  problems  in  th^.  majority  of  cases  are  too  complicated  to 
allow  an  exact  solution  in  the  strict  mathematical  sense.  Therefore,  practicing 
engineers  must  either  replace  the  actual  structure  in  his  analysis  by  a  simpler  one 
which  can  be  solved  or  have  recourse  to  an  approximate  method  of  analysis.  In 


*A  complete  listing  of  references  is  given  on  page  21. 
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recent  years,  numerical  methods,  augmented  by  the  powerful  finite  element  technique, 

2  3 

provide  a  realistic  approach  to  many  complex  thermal  mechanical  problems.  * 

There  is  little  doubt  that  good  finite  element  codes  exist,  but  users  must  be 
cautious  that  the  performance  of  these  codes  is  far  from  uniform.  Generally,  be¬ 
cause  finite  element  codes  are  not  subject  to  regulations,  there  is  no  standard  or 
warranty  on  their  performance.  In  a  recent  report,  a  study  was  conducted  on  the 
thermal  effects  of  an  elementary  type  beam-column  frame  subjected  to  a  temperature 
rise  which  varies  linearly  along  the  beam  and  virtually  across  sections  normal  to 
the  beam.  These  thermal  loads  were  superimposed  with  mechanical  loads.  A  compari¬ 
son  was  made  of  the  results  derived  from  the  application  of  major  finite  element 

4 

systems  (using  beam-type  elements  only)  with  those  from  hand  calculations.  It  is 
not  surprising  to  discover  that  some  well-advertised  programs  could  not  solve  the 
problem  and  others  handled  the  temperature  effects  in  various  ways  which  resulted  in 
a  large  discrepancy  in  the  computed  answers. 

All  major  finite  element  systems,  commonly  known  as  general  purpose  computer 
programs,  possess  certain  thermal  analysis  capabilities.  Most  of  them  have  dem¬ 
onstrated  their  heat  transfer  capabilities  by  providing  solutions  to  radiation, 
convection,  or  heat  conduction  problems.  In  addition,  problems  involving  thermal 
distortion  have  rarely  received  much  properly  documented  attention. 

A  finite  element  code  to  account  for  thermal  expansion  in  a  solid  was  developed 
for  the  in-house  research  project  "Tribology  of  Sliding  Surface  Bearings."  The 
program  is  based  on  a  two-dimensional  model*  using  a  second,  or  higher  order 
interpolation  function  in  the  element  space  that  will  allow  a  diverse  temperature 
field  such  as  a  linear  or  any  nonlinear  temperature  gradient  to  be  prescribed  in  a 
solid  body. 

This  code  has  a  definite  advantage  over  certain  commercially  operated  finite 
element  systems  which  accept  only  constant  or  averaged  temperature  field.  With  the 
new  capabilities,  many  practical  problems  can  be  readily  solved,  for  example,  the 
hot  spot  in  the  landing  deck  due  to  the  concentrated  heat  load,  such  as  those 
generated  by  high  temperature  jet  exhaust,  can  be  more  realistically  represented  by 
the  elements  of  current  development.  Element  mesh  size  and  input  data  set  shall  be 
more  manageable. 


*A  full  3-D  finite  element  model  is  planned  for  the  second  phase  of  the  current 
development . 
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SLIDING  SURFACE  BEARINGS 

Figure  1  illustrates  a  crowned  tilt  pad  bearing"*  situated  under  a  hydrodynamic 
film,  h,  of  lubricant.  The  purpose  of  the  current  study  is  the  development  of  a 
computational  program  using  the  finite  element  technique  to  predict  surface  dis¬ 
tortions  for  pivot-pad  thrust  bearings  as  a  function  of  applied  pressure  and 
temperature  gradient.  Eventually  the  distortion  analysis  will  form  an  integral  part 
of  a  thermoelastic  hydrodynamic  computational  scheme  that  will  have  the  capability 
to  predict  the  behavior  of  lubrication^’^  even  at  its  extreme  conditions,  such  as 
the  boundary  lubrication. 

RUNNER  u 


Figure  1  -  A  Crowned  Slider  Bearing  Pad 
(From  figure  by  S.  Abramovitz) 


The  performance  of  a  hydrodynamic  thrust  bearing  or,  more  specifically,  its 
load  carrying  capacity  (i.e.,  the  sum  of  the  distributed  film  pressure  over  a  bear¬ 
ing  pad),  depends  to  a  large  measure  on  the  profile  of  the  lubrication  film  which, 
in  turn,  is  a  function  of  the  sliding  velocity  vector,  the  thermal  coupling  between 
the  fluid  film  and  the  surrounding  bodies,  the  thermal  and  elastic  distortion  or 
the  crowning  of  the  bearing  surface,  and,  of  course,  the  properties  of  the  lubri¬ 
cant.  The  bearing  sits  tight  on  its  p  ot  against  its  runner,  which  is  a  rotating 
collar  and  operates  satisfactorily  because  it  is  sufficiently  flexible  to  distort 
elastically  and  thermally. ^  The  thermal  distortions  are  of  over-riding  importance 
for  these  bearings. 
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These  assertions  outlined  previously  were  substantiated  by  both  theoretical  and 
experimental  evidence.  It  was  known^  that  the  load  capacity  of  a  truly  flat 
centrally-pivoted  pad  operating  in  air  is,  theoretically,  zero.  The  variation  of 
viscosity  (and  density  as  well)  of  a  liquid  lubricant  with  temperature  accounts  for 
only  a  small  part  of  the  load  capacity  of  such  a  pad  operating  in  a  liquid. 

Convexity,  or  the  surface  crown,  inadvertantly  obtained  in  machining  and  by  load  and 
temperature  deflections  in  service,  probably  accounts  for  a  large  part  of  the 
observed  load  capacity  of  a  centrally-pivoted  pad.  By  controlling  convexity,  the 
load  capacity  of  such  bearings  may  be  increased  appreciably. 

The  standard  procedure,  used  by  lubrication  engineers  to  analyze  the  distortions 

of  the  pad  and  thrust  collar  to  give  film  shape,  is  solving  the  biharmonic  equation 
9  6  8 

of  plate  flexure;  and  the  finite  difference  approach  ’  is  normally  followed.  In 
this  study,  we  set  out  a  scheme  that  uses  versatile  continuum  finite  elements.  It  is 
envisaged  that  the  new  approach  will  remove  all  restrictions  inherited  from  the 
plate  theory,  such  as  shear  deformation  in  a  thick  plate,  etc.  Besides,  the  finite 
element  solution  will  be  advantageous  in  providing  a  precise  representation  of  the 
complex  pad  geometry  and  temperature  distribution,  and  the  solution  is  numerically 
stable . 

In  a  series  of  group  conferences,  the  classical  theory  of  thermoelastic  hydro- 
dynamic  analysis^’ ^  of  the  thrust  bearing  mechanism  was  evaluated,  attention  was 
focused  on  numerical  techniques  for  the  solution  of  finite  bearings,  and  other 
alternative  approaches  were  discussed.  The  major  objective  of  the  joint  research 
program  was  to  develop  a  reliable  tool  to  predict  the  behavior  of  the  thrust  bearings 
in  a  wide  range  of  operating  conditions,  i.e.,  in  the  hydrodynamic  and  boundary 
lubrications . 


THERMOELASTIC  FORMULATION 
THEORY  OF  A  FINITE  ELEMENT  SOLUTION 

When  a  solid  is  under  the  environment  of  a  thermal-mechanical  loading,  the  total 
strains  at  each  point  of  a  heated  body  are  made  up  of  two  parts.  The  first  part  is 
a  uniform  expansion  proportional  to  the  temperature  rise,  0  .  Since  this  expansion 
is  the  same  in  all  directions  for  an  isotropic  body,  only  normal  strains,  and  no 
shearing  strains,  arise  in  this  manner.  If  the  coefficient  of  linear  thermal 
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expansion  is  denoted  by  a,  this  normal  strain  in  any  direction  is  equal  to  a6fi. 
Hence,  initial  strains,  due  to  thermal  expansion  in  the  case  of  a  plane  stress 
model,  become 


thermal  strain  vector 


(1) 


The  second  part  comprises  the  strain  required  to  maintain  the  continuity  of  the 
body  as  well  as  those  arising  because  of  external  loads.  These  strains  are  related 
to  the  stresses  by  means  of  the  usual  Hooke's  law  of  linear  elasticity.  The  total 
strains  are  the  sum  of  these  two  components  and  are  expressed^' ^  as  follows: 


e  =  [a  -va  ]  +  a©  ^ 

x  E  x  y  1 


£ 

y 


1 

E 


[a  -va  ]  +  a0 

y  x 


(2) 


xy 


G  Jxy 


where  G  =  E/2(l+v),  and  E  and  V  are  elastic  constants  commonly  known  as  Young's 
modulus  and  Poisson's  ratio,  respectively.  For  the  range  of  temperatures  encoun¬ 
tered,  e.g.,  bearing  surface  temperature,  the  material  properties  (a  and  E)  are 
taken  as  constant  (temperature  independent)  which  appears  to  be  adequate  for  current 
applications  (see  examples). 

3  12 

Following  a  typical  finite  element  format,  ’  derived  from  energy  consider¬ 
ations,  the  basis  of  computation  for  an  element  is  obtained  in  matrix  notation. 


{F)e  =  f  [B]T[D][B]  dv{6}6  -  {F  }6  -  (F  }e 

J  £  p 


(3) 
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where  { F} e  is  a  vector  of  applied  nodal  forces 

{ F  } e  =  |  [B]T[D]{e  }  dV  is  the  nodal  force  vector  due  to  initial, 

Jy  °  imposed  strain  vector  { eQ } 


(Fp)e  -  J  [N]{p) 
V 


dV  is  the  nodal  force  vector  due  to  distributed 

loads  {p } 


and 


j  IB]T[D] [B]  dV  =  [SK]e  is  the  element  stiffness  matrix 
V 


6  T 

{6}e  =  [u^,v1,u2,V2> . . . ]  is  an  element  displacement  vector 

(ui) 

=  (6 ± } e  where  6±  =  l  i  =  1,  ....  NNPE 

K 


Then,  NNPE  is  the  number  of  nodes  per  element,  for  instance,  NNPE  =  8.  The  inte¬ 
gration  is  normally  carried  out  numerically  as  a  summation  of  quadrature,  such  as 
Legendre-Gauss  quadrature.  A  three-point  integration  along  each  coordinate  line 
appears  adequate. 

Consider  a  general  quadrilateral  element  having  an  arbitrary  shape.  Figure  2 
depicts  an  8-node  element;  its  corner  node  numbers  are  1  to  4  and  the  midside  nodes 


PU.q) 


Y  l.  n  CURVILINEAR 
COORDINATES  LOCAL 
TO  THE  ELEMENT  SPACE 

X 


ELEMENT  ATTRIBUTES 


X 


ELEMENT  EXPANSION 


Figure  2  -  A  Curved  Quadrilateral  Element  for  Thermal  Analysis 
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are  5  to  8.  These  nodes  define  the  boundary  of  the  element  space.  This  element, 

3  12 

commonly  known  as  an  isoparametric  element,  *  is  chosen  for  the  current  application 
for  its  proven  performance  and  simplicity  by  which  the  temperature  field,  TEMP(x,y) 
in  an  element  space,  can  be  described. 


TEMP (x, y)  =  Nx  TEL  +  N2  TE2  +  . . .+  Ng  TEg 
=  TEMP(C.n) 


(4) 


where  TE(I)  is  the  temperature  at  node  I,  and  N  =  N  (£,n)  is  the  interpolation 

3  11 

functions  (or  shape  functions)  where  £  and  n  are  curvilinear  coordinates  local  to 
an  element.  In  this  example,  a  piecewise  linear  or  a  piecewise  quadratic  temperature 
field  can  be  admitted. 

The  material  components  of  a  body  will  expand  or  contract  as  the  temperature  in 
the  neighborhood  of  that  material  point  rises  or  falls  with  respect  to  a  reference 
temperature  TREF,  also  referred  to  as  the  ambient  temperature.  That  is, 

0 (x,y)  =  TEMP (x, y)  -  TREF 


As  the  body  of  a  continuum  is  idealized  into  finite  elements  and  their  connecting 
nodes,  a  set  of  nodal  forces  {F£l  representing  traction  or  compression  between  the 
interface  of  adjoining  elements  will  be  generated  from  Equation  (3). 


{Fe} 


tfx  j  r1  r1 

>=  t  I  [B  ]T[D]{£  }  | J |  d£dn 

TFYi )  j-i  j-i 


(5) 


i  =  1,  NNPE 


These  nodal  forces  may,  at  times,  be  of  a  higher  order  than  those  surface  loads 
which  commonly  occur.  Special  attention  is  called  for  to  ensure  that  these  force 
vectors  are  properly  summed.  Some  detailed  accounts  on  computations  follow. 
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COMMENTS  ON  COMPUTATIONAL  SCHEME 

Consider  the  expansion  load  vector  {F£}e.  In  Equation  (5),  in  which 

3x  3y 

U 

3x  3^ 

3n  Sn 

|j|  is  the  Jacobian  determinant  for  coordinate  transformation.  To  carry  out  this 
calculation,  the  quantities  under  the  integral  of  Equation  (5)  must  be  defined  at 
each  integration  point  •  Begin  with  components  of  strain: 


J  = 


3(?,n) 


T  T 

^xy  3y  3x  3y  —  i  3x  — i 


Now,  we  can  write  the  strain  matrix  [B]  which  relates  the  strain  vector  and  nodal 
displacement  vector  of  an  element 
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(6) 


(e}e  - 

G 

X 

e 

y  1 

.  =  [B]  , 

l“‘| 

3xl 

Jxyj 

3x16 

I”1) 

16xl 

and,  explicitly, 


- 1 

Z 

H* 

3N2 

3N3 

SNa 

0  . 

0  . 

0  ,  . 

.  .  0 

9x  * 

3x  ’ 

9x  ’  ’ 

*•  3x  ’ 

[B]  = 

9N, 

9h’2 

9N3 

3N, 

0  , 

0  , 

0  ,  . 

..  0  , 

1 

1 

3x16 

3y  ’ 

3T’ 

9y  ’  • 

••  9y 

9N. 

3N2 

9N3 

9N8 

9n 

9N2 

9N3 

9N. 

< 

L  W 

w* 

3y"’  ‘ 

“  9y  ’ 

9ir> 

9x~’ 

97“’  * 

••  9x 

The  shape  functions  N ,  =  (C,l)  are  described  in  local  coordinates  to  differentiate 

with  respect  to  global  coordinates,  we  follow  the  standard  rule  of  partial  dif¬ 
ferentiation. 


9N. 

1 


9x 


9N 


i  35 


9N . 


9f,  9x  3n 


3n 

9x 


and 


9N  . 

i 


9y 


3N 

9C 


.  3N  .  , 

i  9C  .  _ i  9n 

9y  9rj  9y 


We  can  proceed  to  form  these  products  by  using  the  Jacobian  matrix  and  its  inverse 
at  each  integration  point 


[Jl 


-1 


in' 

9x  9x 

M  In 

_  3y  9y 
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While  forming  the  Jacobian  determinant,*  note  that 


I  J|  „  9*  h.  _  lx  1* 

|J|  95  9n  9n  9n 


9Na  9NT 
W  3n 


1 


and  that  the  matrix 


9N 

T 

bn! 

— i 

95 

3n 

9Nf  9NT 

9TT  h 


[Qij] 


(7) 


It  is  seen  that 


[Q^j ]  is  skew  symmetric,  i.e., 


In  addition, 


Q(I, J)  =  -  Q(J,I)  and 

Q(I,I)  =0  for  I,  J  *  1.8 


X*  Q(ij) 

|J| 


_  S2  (j ) 


sum  on 


(8) 

(cont. ) 


*To  have  an  inverse  of  the  Jacobian  matrix  implies  |j|  I1  0,  because  the 

12 

Jacobian  index  is  characterized  on  an  element  geometry.  The  user  of  such  an 
element  should  not  allow  the  quadralateral  to  be  in  a  degenerated  form. 
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and 


£  ?(ij> 

i  j  i 


sum  on  1  \ 


SKJ) 

|J| 


(8) 


The  Jacobian  determinant  is  simply 


8 

Z  s2<j>  yi  ■ 

j=l 


u 


Having  the  element  strain  matrix  [B]  defined  term  by  term,  we  can  summarize  the 
expression  for  the  thermal  force  vector  as 


From  Equation  (5), 


ITFX(8) 
TFY (8) 

=  TEMPF  (8,2) 


Sl(l)  0  S2 (1) 
Sl(2)  0  S2 (2) 
SI (3)  0  S3 (3) 


r1  r1  o 

• 

T*1 

C 

o 

_ 1 

fa6  ’ 

(8,2)  =  I  Iii 

SI (8)  0  S2  (8) 

(EV\ 

v  1  0 

< 

a0 

-1  lJl 

0  S2  (1)  Sl(l) 

0  0  “ 

l  oj 

0  S2  (2)  SI (2) 
0  S2 (3)  SI (3) 


0  S2 (8)  SI (8) 


|  J | d£dn 


(9) 
(cont . ) 
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•zz 

tr  iC 


'si(i) 

,S2 (i) 


16x1 


(1+v)  a0  *  DKT  *  H ^ 

i  =  1,8 


(9) 


is  the  weighting  coefficient  corresponding  to  position 
along  Gaussian  quadrature  points  £  ,  » 

respectively,  y 

E  *  t  *  FACT  (NEL) 

(1-v2) 

t  *  FACT  (NEL)  is  the  thickness  of  the  element  number  NEL 
These  thermal  expansion  forces,  TEMPF,  are  the  result  of  temperature  rise,  0, 
inside  an  element  space.  They  must  be  properly  summed  to  form  the  thermal  load 
vector,  TLOAD,  for  the  whole  structure  system,  i.e., 

TLOADX  (NOD)  =  Z  TFX  (NI)  \ 

TLOADY  (NOD)  =  Z  TFY  (NI)  j  SUm  °n  elements 

where  NI  is  a  local  node  number  attached  to  each  individual  element;  for  a  given  NI 

there  is  a  corresponding  global  node  number  NOD.  For  each  node,  the  summation  shall 

extend  to  all  elements  that  are  connected  to  that  node. 

With  these  strain  matrix-elements  completely  described  at  each  integration 

0 

point,  the  element  stiffness  matrix  [SK]  can  be  immediately  computed  in  much  the 

12 

same  fashion.  Finally,  all  those  element  matrices  are  assembled  for  the  complete 

structure.  On  imposing  appropriate  boundary  conditions,  a  system  of  linear  equations 

expressing  the  load  displacement  or  equilibrium  relation  for  the  structure  results. 

12  13 

The  system  is  then  solved  by  the  frontal  solution  technique. 


where  H . ,  H . 

m 

DKT  = 


APPLICATIONS  OF  THE  THERMAL  DISTORTION  CODE 
Preliminary  evaluations  of  the  computed  results  on  the  element  stiffness 
characteristics,  as  well  as  the  nodal  expansion  force  vector,  for  some  basic  element 
configurations  and  temperature  profiles  indicate  that  the  present  program  has  per¬ 
formed  well  in  speed  and  reliability.  In  what  follows,  a  few  selected  examples  are 
given  for  illustration. 

Thermal  Bending  of  a  Cantilever  Beam 
Thermal  Distortion  of  a  Tilt  Bearing  Pad 
Structural  Responses  of  a  Ship  Deck  to  Heat  Load 
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THERMAL  BENDING  OF  A  CANTILEVER  BEAM 


A  simple  cantilever  beam  (Figure  3)  was  subjected  to  a  linearly  varying 

temperature  gradient  across  the  depth  of  the  beam;  no  mechanical  loads  were  applied. 

14 

This  example  is  taken  from  the  ADINA  theoretical  manual  sample  problem  number 

seven  in  which  the  beam  was  modelled  using  three,  16-node  solid  elements.  The 

analysis  was  carried  out  for  material  nonlinearities  only,  and,  by  using  appropriate 

displacement  boundary  conditions,  only  the  portion  of  the  beam  above  the  neutral 

surface  was  included  in  the  finite  element  model.  Partial  results  of  ADINA's 

calculation  are  reproduced  in  Figure  4  showing  the  displacement  response  of  the 

centilever  neutral  surface.  Excellent  agreement  with  the  classical  solution  by 
1 

Boley  and  Weiner  was  obtained. 

The  ADINA's  sample  problem  is  used  here  to  test  the  current  development.  For 
the  present  calculation,  the  complete  beam  is  represented  by  three  second-order 
quadrilateral  elements.  A  constant  thermal  expansion  coefficient,  a  =  0.000006  in./ 
in./°F,  is  used  instead  of  a  slightly  lower  coefficient,  't  =  ot  (TEMP),  that  was 

3  3 

used  in  the  ADINA  example.  As  expected,  the  end  deflection  is  slightly  higher  by 
the  present  calculation,  6^  =  0.0054  in.  (against  0.0052  in.  by  ADINA).  The  complete 
displacement  response  of  the  sample  beam  problem  is  superimposed  on  Figure  4  (de¬ 
flections,  <5  is  measured  along  -Z  direction).  The  performance  of  the  current 
development  appears  favorably  in  terms  of  speed  and  reliability. 

THERMAL  DISTORTION  OF  A  TILT  BEARING  PAD 

The  heat  energy  was  generated  by  the  flowing  lubrication  oil  squeezed  between 
the  sliding  bearing  surfaces.  Part  of  the  heat  was  absorbed  by  the  solid  bearings 
and  the  balance  was  carried  away  through  convection  of  the  lubrication  fluid  into 
the  oil  grooves.  For  the  thrust  bearing  at  given  operating  conditions  of  speed  and 
load,  the  temperature  distribution  within  a  pad  will  converge  to  a  pattern  satisfying 
the  equation  of  heat  balance.  To  demonstrate  the  performance  of  the  finite  element 
code,  a  set  of  input  temperatures,  simulating  a  bearing  pad  under  a  test  situation, 
is  chosen  as  input  data  (Tables  8  and  9  from  Reference  15).* 

The  dimensions  and  material  properties  of  the  bearing  pad,  together  with  a  set 
of  temperature  data  are  shown  in  Figure  5.  Four  8-node  quadrilateral  elements 
connected  to  a  set  of  23  nodal  joints  are  employed  for  its  representation.  The 

*The  pad  geometry  follows  the  description  of  T.  Daugherty's  memo  (Code  2832) 
dated  4  Dec  1980. 
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ENLARGED  VIEW  OF  MESH 
70° F  200° F 


DISPLACEMENT 
CONSTRAINTS  AT 
NEUTRAL  AXIS 


E  30  *  10®  28  x  10® 

v  0.0  0.0 

a  5.5  x  10~®  R0  x  10“® 


lb/in.Z 


in./in./’F 


REFERENCE  TEMPERATURE  -  125°F  -  Tc 


Figure  3  -  Finite  Element  Mesh  of  a  Cantilever  Beam 


Figure  4  -  Elastic  Displacement  Response  of  Cantilever 


184°F,  CONSTANT  ACROSS  THE  TOP  SURFACE 


MATERIAL  PROPERTY:  E  =  29,000  ksi  1/  =  0.30 
TREF  =  70°  F  =  6.0  x  10~6 


Figure  5  -  Bearing  Pad  Geometry  and  its  Finite  Element  Representation 

(With  temperature  input) 

top  surface  is  subjected  first  to  a  constant  temperature  of  184°F  shown  on  line  a, 
and  then  to  a  variable  temperature  increasing  from  166°F  to  184°F  along  the 
direction  of  flowing  lubricant  on  line  b.  For  the  bottom  face,  a  constant  tempera¬ 
ture  of  160°F  is  maintained  in  both  cases.  The  temperature  gradient  is  assumed 
uniform  across  the  thickness  of  the  pad,  otherwise  temperature  varies  linearly  be¬ 
tween  corner  nodes  in  each  element.  The  base  line  of  the  finite  element  model  is 
supported  against  vertical  movement,  but  the  pad  model  is  allowed  to  expand 
laterally  along  the  OX-axis  and  vertically  parallel  to  the  OY-axis. 

Selected  numerical  results  of  surface  distortions,  namely,  displacement  compo¬ 
nents  Y,  are  plotted  in  Figure  6.  Curves  (a)  predict  the  situation  when  the  pad  is 
subjected  to  a  uniform  temperature  gradient  across  the  pad  width.  Here  the  thermal 
distortion  consists  of  a  uniform  thermal  expansion  and  a  simple  thermal  bending. 

The  effect  of  thermal  bending  is  a  symmetrically  crowned  surface  having  a  maximum 

offset  at  the  midspan  of  0.00060  -  0.00045  =  0.00015  in.  which  agrees  very  well  with 

9 

the  hand  calculation  of  0.000155  in.  based  on  the  thermal  bending  formula.  Curves 
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Figure  6  -  Thermal  Distortions  of  a  Bearing  Pad 

(b)  indicate  the  distorted  surfaces  of  the  pad  under  a  skewed  temperature  load,  line 
b  in  Figure  5.  As  a  result,  the  surface  crowning  shifts  towards  the  trailing  edge. 
For  the  top  surface,  a  crown  height  of  0.00010  in.  occurs  about  1.41  in.  from  the 
leading  edge. 

STRUCTURAL  RESPONSES  OF  A  SHIP  DECK  TO  A  HEAT  LOAD* 

A  ship's  deck  can  vary  in  size  and  shape,  and  differ  in  the  method  of  con¬ 
struction  and  design  arrangement.  However,  typically  the  deck  consists  of  a  welded 
steel**  plate  roughly  1/2-  to  1-in.  thick  depending  on  its  functional  requirement. 
The  plate  is  cross  stiffened  longitudinally  and  transversely  by  Tees  ranging  in  size 
from  10  to  18  in.;  these  T-beams  are,  in  turn,  framed  into  stiffened  wall  panels. 
Further,  a  system  of  bulkheads,  stanchions,  and  intermediate  wall  panels  were 
installed  to  augment  the  load-carrying  capacity  of  a  deck  structure.  It  is  clear 
that  the  combination  of  these  support  systems  results  in  deck  construction  which  is 


*This  problem  was  funded  by  the  surface  ship  structure  block  program.  Task  Area 
SF  43  400  391  and  Work  Unit  1730-035. 

**Alternatively  aluminum  constructions  were  employed. 
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highly  redundant.  Presently,  our  interest  is  to  identify  the  primary  structural 

behavior  of  a  deck  subjected  to  a  thermal  loading.  We  shall  focus  our  attention  on 

a  basic  transverse  beam-column  structure  unit. 

Figure  7  shows  a  transverse  bent  composed  of  a  3/4-in.  thick  deck  plate 

stiffened  by  an  18  in.  by  0.5  in.  web  plate  forming  a  Tee  section  girder  having  a 

4 

moment  of  inertia,  I  =  488  in.  The  column  has  a  cross  section  of  12  in.  by  0.5  in. 

®  4 

of  which  the  moment  of  inertia,  I  =  72  in.  The  bent  has  an  outside-to-outside 

’  c 

width  of  28  ft  and  head-room  of  8  ft  to  the  deck  below. 

For  simplicity,  we  set  up  symmetric  temperature  data  which  somewhat  resemble 
the  heating  load  impinging  upon  a  HELO-DECK  due  to  the  exhaust  gas  from  a  moderate¬ 
sized  V/STOL  (Vertical,  Short  Take  Off  and  Landing  Craft)  such  as  the  AV8  Harrier. 
Temperature  could  rise  to  a  much  higher  level  if  the  deck  is  exposed  to  the  jet 
exhaust  of  a  high  performance  major  V/STOL  craft. 

A  maximum  surface  temperature  of  250°F  occurs  at  midsection  near  the  heating 
source.  It  drops  off  quickly  to  130°F  at  a  distance  35  in.  to  50  in.  away;  then 
gradually  reduces  to  95°F  near  the  edge  of  the  deck.  The  under  side  of  the  deck 
plate  has  a  maximum  temperature  of  170°F  at  the  midspan  and  tapers  off  gradually  to 
90°F  near  the  edge.  Approximately  linear  temperature  variation  is  assumed  across 
the  thickness  of  the  deck  plate.  Other  parts  of  the  structural  components  will  have 
an  ambient  temperature  of  70°F. 

For  a  symmetric  heating,  we  need  to  consider  only  half  of  the  transverse  bent 

which  can  be  represented  by  twelve  8-node  elements*  and  a  total  of  56  nodal  joints 

in  Figure  8.  Let  us  assume  the  midsection  of  the  bent  rests  on  top  of  a  bulkhead 

which  provides  a  solid  support.  If  we  leave  the  bottom  of  the  bent  free  to  move,  it 

will  have  a  distorted  shape,  as  shown  in  Figure  9.  At  node  54,  the  center  of  the 

54  54 

bottom  edge,  there  is  a  movement  of  0.4628  in.  (AX  =  0.3143  in.,  AY  =  0.3391 
in.)  and  a  clockwise  rotation  of  0.0031  radians.  If  we  move  the  same  point,  node 
54,  back  to  its  original  position,  a  reactive  force  of  3.422  kip  (Rx  =  -0.106  , 

k 

Ry  =  3.42  )  is  required  to  maintain  the  structural  continuity,  i.e.,  between  the 
bent  structure  and  its  supporting  deck  (01  level). 


*A  more  elaborate  finite  element  model  composed  of  120  elements  or  more  can  be 
set  up  when  the  loading  situation  calls  for  such  an  analysis. 


95°  F 


DEFORMED  SHAPES  REPRESENTED  BY  A  NEUTRAL  SURFACE 


o™4  P 


sj7  =  0.0186 


3.42' 


Figure  9a  -  Base  is  Free 


Figure  9b  -  Base  is  Supported 


Figure  9  -  Distortions  of  a  Half  Transverse  Bent  by  Deck  Heating 


CONCLUSIONS 

A  thermal  distortion  analytical  capability  has  been  implemented  in  a  finite 
element  computer  code.  The  numerical  results  have  correctly  predicted  the  structural 
responses  to  a  variety  of  thermal  loadings.  It  is  known  that  thermal  distortions 
are  of  over-riding  importance  to  the  performance^  of  a  tilting  pad  bearing.  The 
present  code  capable  of  representing  precisely  both  a  linear  and  even  an  arbitrary 
temperature  distribution,  furnishes  an  imperative  feature  to  the  execution  of  a 
joint  task  on  the  thermoelastic-hydrodynamic  calculation  of  sliding  surface  bearings. 

Recent  studies  on  thermal  structural  responses  also  indicate  some  inconsis- 

4 

tencies  in  the  thermal-structural  capabilities  of  some  major  finite  element  systems 
that  have  been  available  commercially.  The  present  development,  even  though  set  up 
on  a  linear  thermoelastic  framework  with  temperature  independent  material  properties, 
can  be  adopted  to  identify  shortcomings  that  may  exist  in  certain  major  finite 
element  systems.  In  addition,  the  present  code  will  provide  an  independent  means 
capable  of  assessing  rapidly  a  large  number  of  structural  responses  under  the 
prescribed  thermal  environments  to  which  only  a  numerical  solution  offers  a  viable 
approach. 
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11 

W. 

Ellsworth 

1 

15 

W. 

Morgan 

1 

154 

J. 

McCarthy 

2 

1542 

J. 

Bai 

1 

17 

w. 

Murray 

1 

1706 

c. 

Eynon  (m) 

1 

173 

A. 

Stavovy 

1 

1730.5 

J. 

Adame hak 

27 

1730.5 

J. 

Ma 

1 

174 

I. 

Hansen 

1 

177.4 

M. 

Wright  (m) 

CENTER  DISTRIBUTION  (continued) 


Copies 

Code 

Name 

1 

1844 

S.  Dhir 

1 

28 

J.  Belt 

1 

283 

1 

2832 

J.  Dray 

5 

2832 

T.  Daugherty 

10 

5211.1 

Reports  Distribution 

1 

522.1 

Library  (C) 

1 

522.2 

Library  (A) 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM 
INAHY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MtMuH ANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY-CASE 
BASIS. 


